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We investigate non-Gaussianity in general multiple-field inflation using the formalism we developed 
in earlier papers. We use a perturbative expansion of the non-linear equations to calculate the 
three-point correlator of the curvature perturbation analytically. We derive a general expression 
that involves only a time integral over background and linear perturbation quantities. We work out 
this expression explicitly for the two-field slow-roll case, and find that non-Gaussianity can be orders 
of magnitude larger than in the single-field case. In particular, the bispectrum divided by the square 
of the power spectrum can easily be of 0(1-10), depending on the model. Our result also shows 
the explicit momentum dependence of the bispectrum. This conclusion of large non-Gaussianity is 
confirmed in a semi-analytic slow-roll investigation of a simple quadratic two-field model. 

I. INTRODUCTION 

The assumption that inflationary fluctuations are Gaussian is a good starting point for the study of cosmological 
perturbations, but it is only true to linear order in perturbation theory. Since gravity is inherently non-linear, and 
most inflation models have (self-)interacting potentials, non-linearity must be present at some level in all inflation 
models. Hence the issue is not whether inflation is non-Gaussian, but how large the non-Gaussianity is. With 
increasingly precise CMB data becoming available in the near future from the WMAP and Planck satellites and other 
experiments, we might well hope to detect this non-Gaussianity. This would offer us another key observable to help 
constrain or confirm specific inflation models and the underlying high-energy theories from which they are derived. 
As a rough order of magnitude estimate, we note that non-Gaussianity will be detectable by Planck if the bispectrum 
(the Fourier transform of the three-point correlator) is of the order of the square of the power spectrum . 

It follows that to compute the predicted amount of non-Gaussianity in specific inflation models we need to go 
beyond linear-order perturbation theory. In 0.0] we introduced a new formalism to deal with the non-linearity during 
inflation. We will not again summarise the other work dealing with this subject, references for which can be found in 
or a recent review Q . Our formalism is distinguished by being based on a system of fully non-linear equations for 
long wavelengths, while stochastic sources take into account the continuous influx of short-wavelength fluctuations 
into the long- wavelength system as the inflationary comoving horizon shrinks. The variables used incorporate both 
scalar metric and matter perturbations self-consistently and they are invariant under changes of time slicing in the 
long- wavelength limit. 

The advantages of our method are threefold: (i) it is physically intuitive and relatively simple to use for quantitative 
analytic and semi-analytic calculations; (ii) it is valid in a very general multiple-field inflation setting, which includes 
the possibility of a non-trivial field metric; and (Hi) it is well-suited for direct numerical implementation. The first 
point was already demonstrated in where we computed the non-Gaussianity in single-field slow- roll inflation, while 
the third point is the subject of a forthcoming paper [(|. The present paper is dedicated to the exploration of the 
second point, as well as the first. 

In [I| we found, confirming what was known in the literature beforehand (see e.g. Q), that non-Gaussianity in the 
single-field case is too small to be realistically observable, because it is suppressed by slow-roll factors (actually the 
scalar spectral index n — 1). However, there have been long-standing claims in the literature (see e.g. 0, 0) that 
specific multiple-field models can, in principle, create significant non-Gaussianity. Indeed, there has been growing 
recent interest in models which can produce large primordial non-Gaussianity [T(il ] . A feature shared by most of these 
models though is that this non-Gaussianity involves some mechanism operating after inflation, usually (p)rehcating 
or later domination of a curvaton field. In this paper we investigate for the first time general multiple-field inflation, 
not just specific models, presenting a method by which to accurately calculate the resulting three-point correlator 
during inflation. We find that it is possible to get significant primordial non-Gaussianity without invoking some 
post-inflationary mechanism even for the simplest two-field quadratic potential. 

The key mechanism for the production of this large non-Gaussianity is the superhorizon influence of isocurvature 
perturbations on the adiabatic mode. The former feed into the latter when the background follows a curved trajectory 
in field space. Note that the example studied in section IV Dl illustrates that there is no need for the potential to be 
interacting. Our aim is to push forward towards a tractable non-Gaussian methodology for the new era of precision 
cosmology which confronts us. 
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This work is organised as follows. In section[H]we give the equations from Q that are used as the starting point for 
the present investigations. In section UTTl we then derive the general solution for the relevant quantities in multiple- 
field inflation, culminating in a general expression for the three-point correlator of the adiabatic component of the 
curvature perturbation, without any slow-roll approximations. This integral solution - equation H29|) - is a very useful 
calculational tool because it gives the three-point correlator entirely in terms of background quantities and linear 
perturbation quantities at horizon crossing. In section llVl we make a leading-order slow-roll approximation to work 
out the various contributions in the general solution more explicitly. Finally in section we calculate the bispectrum 
in an analytic treatment of the two-field case with constant slow-roll parameters. We find that the result can be orders 
of magnitude larger than for single-field inflation. This result is confirmed with a semi-analytic slow-roll calculation of 
an explicit model with a quadratic potential in section IV D\ Our method yields the full momentum dependence, not 
just an overall magnitude, and we find that there can be a difference of the order of a few between opposite extreme 
momentum limits. We conclude in section fvTI Parts of this paper, in particular section |3 are rather technical, so 
some readers might be interested in referring to [TT| | first, which contains a simplified derivation of only the dominant 
non-Gaussian contributions in multiple-field inflation, along with a summarised discussion. 



II. MULTIPLE-FIELD SETUP 

Since in this paper we are explicitly working out the general non-linear multiple-field formalism of , we refer the 
reader to that paper for derivations and more details of the initial equations. Here we just briefly describe the context 
and give the relevant equations and definitions to be used as starting point for further calculations. 

We start from a completely general inflation model, with an arbitrary number of scalar fields (f> A (where A labels 
the different fields) and a potential V(4> A ) with arbitrary interactions. We also allow for the possibility of a non-trivial 
field manifold with field metric Gab ■ We will consider only scalar modes and make the long- wavelength approximation 
(i.e. consider only wavelengths larger than the Hubble radius l/(aH), where second-order spatial gradients can be 
neglected compared to time derivatives) 1 . The spacetime metric and matter Lagrangean are given by 

ds 2 - -N 2 (t, x)dt 2 + a 2 (t, x)dx 2 , C m = —g* 'd^ A G AB d v ^> B - V(<f> A ), (1) 

with a the local scale factor and N the lapse function. The local expansion rate is defined as H = a/(Na), where the 
dot denotes a derivative with respect to t. The proper field velocity is 77 A = <j) A /N, with length 77. We also define 
local slow-roll parameters as 

~u \ k2jj2 ~au n 3Hn A + G AB d B v ~ A . . v B d A vn B ^n A o _ A 

where k 2 = 8irG — 8tt /rn 2 ^ and is a covariant derivative with respect to the field 4> B . For the first part of the 
paper we will not make a slow-roll approximation, and consider these definitions as just a short-hand notation. When 
we do make this approximation, from section llVl e and rj A are first order in slow roll, while £ A is second order. Finally, 
we choose the gauge where 

t = bx(aH) & NH = (1-e) -1 . (3) 

In this gauge horizon exit of a mode, k = aH, occurs simultaneously for all spatial points and calculations are simpler. 

We will make use of a preferred basis in field space, defined as follows. The first basis vector e A is the direction 
of the field velocity. Next, e A is defined as the direction of that part of the field acceleration that is perpendicular 
to e A . One continues this orthogonalisation process with higher derivatives until a complete basis is found. Explicit 
expressions can be found in [jj], here we only give e A = II A /II. Now one can take components of vectors in this basis 
and we define, for example for Q A (defined below in ©) and fj A : 

C = e mA C A , V ll =efvA, t = e A U- (4) 



1 Formally this corresponds with taking only the leading-order terms in the gradient expansion. We expect higher-order terms to be 
subdominant on long wavelengths during inflation, but this statement has only been rigorously verified at the linear level. A calculation 
to higher order in spatial gradients, or, even better, a full proof of convergence of the expansion, would be desirable. See |l2j for more 
details on the validity of the gradient expansion beyond linear theory. 



3 



Note that, unlike for the index A, there is no difference between upper and lower indices for the m. By construction 
there are no other components of fj A , so that one can write fj- 1 — \f] A — f^e A \. We also define 

1 

e m A?>te A , (5) 



NH 



where T> t is the covariant time derivative containing the connection of the field manifold. Z mn is antisymmetric and 
only non-zero just above and below the diagonal, and first order in slow roll. Its explicit form in terms of slow-roll 
parameters can be found in [l3j : here we only need that Z 12 = — Z 2 i = —rj- 1 . 

As discussed in 0, Q it is useful to work with the following combination of spatial gradients to describe the fully 
non-linear inhomogeneities: 

( A (t, x) = e A (t, x)di lna(t, x) - di<j> A (t, x) C = S ml d t In a - -JL e mA d t <j> A , (6) 

y2e(t, a:) y2e 

which is invariant under changes of time slicing, up to second-order spatial gradients 0, IT^ |. Note that, when 
linearised, Q (the m=l component of C™ 1 ) is the spatial gradient of the well-known ( from the literature, the curvature 
perturbation. In [j| we derived a fully non-linear equation of motion for £™ without any slow-roll approximations: 

C - of = s™ 

(3 - 2e + 2jjU - 3e 2 - 4S?" )J mn , 2Z mn ^ Qn , S mn = ^ m (7) 



(1 - e) 2 1 - e y 2 (1 - e) 2 

where 0™ is the velocity corresponding with and 

1 2£ R m iin + {l-e)Z mn + Z mp Z pn +(3-2i+2^-~e 2 -2^)^ (8) 



H 1 n £ V / 1 — e 

+ (3e + 37)11 + 2f~ 2 + 4e~ryll + (t^) 2 + fU) <5 mn - 2? ((3 + e + 2rjU)* ml J n i + ?? X (5 m i<5„ 2 + Wni)) , 

where V mn = e A l (T>BdAV)e R and R m iin = ^m^ABCD^i e \^n with R A bcd the curvature tensor of the field manifold. 
Although for the first part of the paper we will not make a slow-roll approximation, we give here immediately the 
leading-order slow-roll approximation of J7J, which we will be using in the second part, to show that things simplify 
considerably in that case: 

£771 QTrl 



3C + + 3Z mn + 3(e + jjH)5 m „ - 6eS ml S nl j Q = JT< 



The stochastic source terms <S>™ and J™ are given by 
k t d 3 k 

Sr = -^Jf-j j^W(k)Q^ n (k)a n (k)ik ie ^- 



d 3 fc 



o lin mA - n lin (h\ i 1 + e + ¥)$n P + Z np 



1 - e 



a p (A;)ifciC Ifex + c.c., (10) 



where c.c. denotes the complex conjugate. The perturbation quantity is the solution from linear theory for 

the multiple-field generalisation of the Sasaki-Mukhanov variable Q = — av / 2lC/ K - It can be computed exactly 
numerically, or analytically within the slow-roll approximation 13]. The a m {k) are a set of Gaussian complex random 
numbers satisfying 

(a m (k)a* n {k')) = 5 mn d 3 {k~k'), (a m (k)a n (k')) = 0. (11) 

The quantity W(fc) is the Fourier transform of an appropriate smoothing window function which cuts off modes with 
wavelengths smaller than the Hubble radius; we choose it to be a Gaussian with smoothing length R = c/(aH) = ce - *, 
where c « 3-5: 

W(k) = e - k2R2 / 2 W(k) = k 2 R 2 c- k2R2 / 2 . (12) 

Since £™ and 6™ are smoothed long-wavelength variables, the appropriate initial conditions are that they should be 
zero at early times when all the modes are sub-horizon. Hence, 

lim C = 0, lim df = 0. (13) 

t— *■— oo t— > — oo 
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A key aspect of the system or © is that it is fully non-linear. All functions in the coefficients on the left-hand 
side of the equation, like e(t, x), and in the sources on the right-hand side are inhomogeneous and depend on C,™ and 
0™ via a basic set of three constraint equations: 



di In a = 


-dilnH 


d4 A = 




v t n A = 


-- — H 

K 



1-6 

€ 



1 - e 



1 - e 



(14) 
(15) 
(16) 



Using only these three constraints one can compute the spatial derivative of all relevant quantities, keeping in mind 
that Q? = e mA 6f - (1 - t)~ x Z mn Q. Note that in our gauge W depends on t only and does not get any non-linear 



contributions. 



III. GENERAL ANALYTIC SOLUTION 



In this section we investigate how to solve the system Q analytically and give formal expressions for the solution. In 
the next sections we will investigate cases where we can determine the solution more explicitly. We start by rewriting 
the system (jTJ) into a single vector equation: 



Via(t,X) + Aab(t,x)v ib (t,x) = b ia (t,x), 



'P 



J? 

V ; J 



(17) 



Here the indices a,b, . . . label the components within this 2A-dimensional space (with N the number of fields). The 
matrix A can be read off from Q and has the following form: A2 m -i,2m — — 1, ^2m,2n = B mn and A2 m ,2n-i = 
2 mn /(l — e) 2 , where Q mn is the matrix between parentheses in the second equation of Q and m,n — 1, 2, . . . , N . All 
other entries of A are zero. 

Equation (|17l) is non-linear since the matrix A(t, x) and the the source term bi(t, x) are inhomogeneous functions 
in space and depend on the Vi through (|X4|> — (|lfc>|) . It can be solved perturbatively as an infinite hierarchy of linear 
equations with known source terms at each order (see also |j). We expand the relevant quantities as 



(1) , (2) 



Then the equation of motion for is 



b l = b^+b^ + 

(m) 



A(t, x) = A (0) (t) + A {1) (t, x) + A (2) (t, x) + 



(18) 



.(to) 



'(t,x), 



lim 



( m ) n 
V ia = 0, 



lim) _ Am) 



4("»-2) 0') 



(19) 



Let us recapitulate the meaning of the various indices, to avoid confusion. The index i = 1,2,3 labels the components 
of spatial vectors. The indices A,B,... — 1, . . . , N label components in field space. These indices will not occur in 
the rest of the paper, since they have been replaced by the indices m, n, . . . = 1, . . . , N that label components in field 
space within the special basis as defined in (0J. Next, the indices a, b, . . . = 1, . . . , 2A label components within the 
2A-dimensional space consisting of both £ and 9 as defined in (|17|) . Finally there are the labels within parentheses 
that denote the order in the perturbative expansion defined above. Only with the i and A, B, , . . is there a difference 
between upper and lower indices. 

We now show that the source term b[ m ^ is known from the solutions for «j up to order (m— 1). The equation for i;^ 
is linear by construction: A^ depends only on the homogeneous background quantities, and the only x dependence in 
b^ is in the e lfc x , for the rest it depends on homogeneous background quantities. All of these are in the end functions 
of just H, (j> , and H A via their definitions. To go beyond linear order all these background quantities are perturbed 
as follows (C stands for any of the quantities to be perturbed, for example e, rj- 1 , etc.): 



C(t,x) = C^(t) + C (1) (t,x) + ... = C (0) +d~ 2 d i (d i C) (1) + 



(20) 
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where we use (H3J-(Uni to compute <9;C and is some homogeneous (space-independent) vector that is the result 
of that calculation. Next, to compute one simply repeats this process with the vector C, and continues in this 
way order by order (of course there is also a d~ 2 d l vf} term at second order, etc.). Then it is easy to see that 
b\ m ^ depends only on . . . , uj m x \ and hence is a known quantity at each order. 
The solution of equation 1(19(1 for vj- can be written as 

v£\t,x) = f dt'G ab (t,t')b^(t',x), (21) 

J — oo 

with the Green's function G ab (t,t') satisfying 2 

d 



dt 



G ab (t,t') + A^(t)G cb (t,t r ) = 0, KmG ab (t,t') = 6 ab . (22) 



It is important to note that this Green's function is homogeneous, a solution of the background equation. It has to 
be computed only once, and can then be used to calculate the solution at each order using the different source terms 
as in (|2U. We write explicitly for the first two orders: 

bf}{t,x) = (d-Wvg&xj) J ^W(M)lll(MK(^)*e itl + c.c. (23) 
Comparison with 110(1 shows that X am is given by the following equations: 



y — D lin y — - 

yY 2n — l,m — r— nrrn ^2n,m — 



av Ze ay Ze 



O lul - O 



up 



1 - e 



(24) 



The quantity X amc is derived from X am as in (|20|) . but in addition it also contains the perturbation of the basis 
vector e m inside the a m . In the same way we define A^ b \t, x) — A^ b ' c (t)d~ 2 d l v^ c '(t, x). 

Using the solution |J2U, valid at each order, and the relations 1)11(1 to compute the averages, it is now straightforward 
to write down the general expressions for the two-point and three-point correlators of the adiabatic (m = 1) component 
of £ m = <9 -2 9 l (™, which is the a = 1 component of Vi a , or rather their Fourier transforms, the power spectrum and 
the bispectrum. Making use of the short-hand notation 

dt' G ab (t,t')W(k,t')X^(k,t'), (25) 



we find for the power spectrum: 

(\C^\k,t)\ 2 )=v[]l(kM 1 l*(k,t)+c.c. (26) 

We emphasise here that the quantity vhm defined by 1(25(1 is simply made up of the linear C^mm ^mn mode functions. 
One needs to evaluate the Green's function G ab (t,t') and to perform the integral ((25|) when the linear source terms 

Qnim Qmn m ^fcm are known only up to horizon crossing, k aH. However, where analytic solutions are available 
for k <C aH , or in a fully numerical approach, we can dispense with the integral ((25(1 by using the linear solution on 
super-horizon scales. 

To find the bispectrum the calculation is slightly longer. One first has to compute 1 (t , x) . As we noted in 
the single-field case in [5|, (C ) is indeterminate. To remove this ambiguity and also require that perturbations 
have a zero average, we define ( m = ( m — (0™)- Expanding ( m = <^( 1 ) m + Q 2 ) m i the three-point correlator will be 
a combination of the different permutations of (C' 2 ' 1 {t, x\)Q^ 1 (t, x^)^ 1 ' 1 (t, afc)), and the bispectrum is its Fourier 



2 To be precise, the Green's function is actually defined as the solution of 1221 with 8(t — f') on the right-hand side instead of zero. The 
solution is then a step function times what we call the Green's function. This step function has been taken into account by changing 
the upper limit of the integral in 1211 from oo to t. 
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transform. The intermediate steps are given in more detail in the explicit calculation in section[V] here we go directly 
to the end result for the bispectrum of the adiabatic component: 



(2) 



C 1 ^, x l )( l {t, a&K'ft aa)> (fci, fc, h) = (2vr) 3 5 3 (fc 1 + k 2 + fc 3 ) f(k u k 2 ) + f(k u k 3 ) + f(k 2 , fc 3 ) 



(27) 



with 



(*', + < J (*', *(*', + ce. + fc - k'. (28) 



If Qm n is real, this simplifies to 



fc+fc'l 2 

t 



x / dt'G la (t,t') i(2 c (M')W(M')-€(«')^(M') «W(fc',i') + ^fc' 



(29) 



This integral expression is a key result of this paper. Using our methodology, the three-point correlator with full 
momentum dependence has been expressed as a single time integral over quantities determined by the background 



and 



,(!) 



model and the linear perturbations, that is, respectively the matrix A ab ' c and the solution Qllf n embedded in X a Jic <|24|) 
5j (in both of which the background is also implicit). Of course, one also has to find the Green's function 



G a b from (|22|l . but, like the equation for Q^ nB in Q, this is a linear ordinary differential equation for which there is no 
serious impediment to finding a numerical solution, in cases where an analytic or semi-analytic solution is unknown. 
The integral solution (|29|) . then, demonstrates that the calculation of the three-point correlator is straightforward and 
tractable. It is, in principle, similar to calculations of the power spectrum, where accurate estimates can be found 
from background quantities, for example, in the slow-roll approximation. Here, we only have to supplement this with 
the amplitudes of the linear perturbation mode functions Q l ™ n (k, t) and the closely related Green's function G ao (t, t'). 
In section IV Dl using this methodology, we provide some quantitative semi-analytic results for the bispectrum of a 
two-field inflation model with a quadratic potential. 

Before closing this section a final comment is in order. A feature of i|29|) which may at first sight cast doubt on 
its utility for quantitative calculations is its apparent dependence on the ad hoc choice for the functional form of the 
window function W(fc). Closer inspection reveals that the second term of (|29(l (the A term) does not depend on W(k) 
for scales sufficiently larger than the horizon. This is evident from the fact that 1)25(1 is simply the solution to linear 
theory smoothed on scales larger than the horizon. Any properly normalised window function with W(fc) — > I for 
scales sufficiently larger than the horizon will produce the same final answer. The A term represents the non-linear 
evolution on superhorizon scales and, as we show below, it describes an integrated effect which can lead to large 
non-Gaussianity. In contrast, the X term arises from perturbations around horizon crossing and may depend on W. 
We find below that this term is localized around horizon crossing and that it does not give rise to observationally 
interesting effects. Section IV El provides more discussion on these points. In ^lj we explicitly show that taking a step 
function instead of a Gaussian as window function does not change the leading-order integrated effects. 



IV. SLOW-ROLL APPROXIMATION 



The perturbation quantity Q^ n can be computed exactly numerically, or analytically within the slow-roll ap- 
proximation where all slow-roll parameters are assumed to be smaller than unity. The latter was done in p!3 | to 
next-to-leading order in slow roll: 3 



2Vk 



e(-) 

\kRJ 



l+D 



(30) 



3 Compared with the solution in Il3l there is an extra factor of 1/ V2- It has to be introduced to take into account the difference between 
the classical Gaussian random numbers a, which have (aa*) = (a* a), and the quantum creation/annihilation operators a T , a, which 
have (at a) = 0. In 5] we introduced this factor of l/\/2 in the analogue of equation 1101 . which leads of course to identical results. 
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where the matrices D and E are defined by 



= eS„ 



2e6 m i5 n i 



Vm 



■AH 2 ' 



E mn = (1 - e)5 mn + (2 - 7 _e - In 2) D r . 



(31) 



with 7s Euler's constant. Overall unitary factors that are physically irrelevant have been omitted. Using this 
expression the source terms are given by 



J? 



K 



d 3 k 

aV2iJ (2tt) 3 /2 

K 



aV2e 



D„ 



VV(fe) Q^ n {k)a n {k) ih e ik x + c.c. 

d 3 A; 



(2i + ^)6 



(2tt) 3 / 2 



Q%(k)a p (k)ik i e 1 ** + c.. 



(32) 



Now when computing X^mc as defined in l|23|) . or any higher-order terms in the perturbative expansion, we would 
in principle have to make the background quantities in Q y ™ n dependent on x and perturb them according to l|2(J|l . 
However, from (|30fl we see that Q l ™ n depends on x only beyond leading order in slow roll (to leading order it is just 



given by Q 



lin 



z/(2k 3 / 2 R)S mn ). Hence in a leading-order slow-roll approximation the only non-linear parts in the 



source terms are the factors in front of the integrals in i|32|l . plus the basis vector e p inside the a p . 

Within the leading-order slow-roll approximation, we now look at the two-field case, to make things a bit more 
explicit. In that case the matrix A in (|17|) is given by 





-i 











3 


— 6r] ± 
















Vo 





3x 


3 



V22 

AH 2 



+ e + ?f 



(33) 



The quantity x is nrs t order in slow roll. Here we used © and the relation, valid to leading order in slow roll (see 



1 ne qui 

3 . g . m 



AH 2 AH 2 



= (e - 77")dml - 77 On 



(34) 



Using the constraints (|14|I -H16 |I we can compute the spatial derivatives that are needed to calculate A. Some of these 
were given in a general form in [3j; to first order in slow roll for the 9 coefficients and to second order for the ( 
coefficients we find in the two-field case, 

(cUna) (1) - -(cUntf)« = -e£ (1)1 , 

(ftlng)W = -2^ 1 -2( e - + ^)c} 1)1 +2^c| 1)2 , (35) 
= ^ 1 + (3 - 3e + ^ 2 - ((e - 2^ + ^) Cf > 1 + (s* + (E + V l] W ] (v^f) Cp, 
(d lX ) {1) = (3 - 5e + 77II )0P 1 - S^- 1 -^ 13 2 - 1 - (rh + (6 - 10e + 2^ + 3 X )tj x ) C 4 (1) 

introducing the two second-order slow-roll quantities tpi and ^2 as short-hand notation: 

Vr ss 2e X +(e-^ ll +3(»r L ) 2 +£ l + — ^5 = 2e ~( x + 2 ? 7 ll )+4(^) 2 -^-i ? (U 11 i-U 221 ), 



i'2 



(lie + 2^H -3x)j? X +f 



k 3iJ 2 

^1^222 



K 3iJ 2 



K 3iJ 2 ' 



(36) 



with U„, 



e^efe^Pc^fl^l''. The reason for this specific definition of t/> 2 will become clear later on. Since we 



have only two fields, the notation £ _L is unambiguous. To compute diX we used that in the two-field case, because of 
the orthonormality of the basis vectors, Vie^ = — e^(e 2 B 2?ief ). All slow-roll parameters in these expressions take 
their homogeneous background values. From this we find that the rank-3 matrix A a b c (defined below (|24l0 is 



A = 



/0 0\ 

-6 f-(e- 2f}W)r)- L - £ ± ,fj- L ,3x + (e + V ll )fj 11 - (t?- 1 ) 2 ^- Al + if^j 



^0 3(-^i,3-5e + r)",-i/' 2 -(6-10e + 27?ll+3x)?y ± ,-3^ ± ) oj 



(37) 
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In the same way we find that the matrices X am and X amc , defined in Q24[l. are given by 



X 

X'22 
X42 




X 



aV2e 



/ (2e + f/H , 1, -fj^, 0) -(ry- 1 , 0, e + 77!!, 1) \ 

277- L (?7- L ,0,e + f?ll,l) X 22 

(^,0,? + ryll,l) (2e + r7U,l,-^,0) 

V-x(^,0,e~ + r7H,l) X 42 / 



2fc 3 / 2 ' 



(38) 



= 2 



(e + Zfj^fj 1 - - t, 2fj x , 3 X + (e + r? 11 )^ 11 - 2(f, ± ) 2 , 3 - 3e + 77U 
(Vi - (2e + ?7 ll )x,-3 + 5e-?7 11 - X , i>i + (6 - 10? + 277 11 + 4x)?? ± ,3rr 



where we used J3DJ|, ^2), and ||3SJ. Note that -(<9;ai)/a 2 = (<9,a 2 )/ai = r^Ci + (e + ^")Cf + to leading order in 
slow roll. 

In general the Green's function G a b(t,t') cannot be expressed in closed form, since the time-dependent matrix A 
does not commute at different times. It can be formally represented as 



G(t,t') =Texp 
where T denotes a time-ordered exponential: 



-J A(s)di 



Texp 



/>t />t pt rS pt />S pS 

/ A(s)ds = 1 - / A(s)ds + ds ds'A(s)A(s') - / ds / ds' / 
Jt' Jt' Jt' Jt' Jt' Jt' Jt' 



ds"A(s)A(s')A{s") + 



(39) 



(40) 



This formal expression is standard in quantum mechanics and quantum field theory (see e.g. |l5j ^l where, viewed as 
a perturbative expansion, the first few terms in the series are kept when the operator A contains a small parameter. 
In our case, however, not all elements of A are first order in slow roll, so that a truncation at any finite order is a 
bad approximation. Moreover, even if A were first order in slow roll, one should still be careful, because the time 
interval in the integrations can easily be of the order of an inverse slow-roll parameter. The time-ordered exponential 
can be written as an ordinary exponential plus terms which contain (nested) commutators. For example, the second 
and third order terms in the series (|40|l can be written as 



t ds£ds'A(s)A(s') = lQ\ 



A(s)ds ) +- ds ds' [A(s), A{s')] 9(s - s') 



(41) 



and 



ds' / ds"A{s)A(s')A(s") = 



3! 



A(s)ds) +-/ ds / ds' / ds"A(s)[A(s'),A(s")]Q(s' - s") 



+ 1 [ ds [ ds' / ds"([A(s),A(s>)\ A(s") - A(s") [A(s),A(s')])q(s - s')9(s - s") (42) 



ds 



respectively, where is the step function, and similarly for higher orders (see |16| for general expressions). 

There are basically three ways to proceed with this expression. In the first place we can, if we are interested only 
in relatively short time intervals, neglect the commutator terms in the expansion of the time-ordered exponential 
and write it as an ordinary exponential. The commutator terms all contain a difference of slow-roll parameters at 
different times, as opposed to the terms of the ordinary exponential that have just a slow-roll parameter at one time. 
Hence, for small time intervals, the commutator terms are a slow-roll order of magnitude smaller. Then we have an 
exact analytic solution in closed form for the Green's function. Secondly, as will be the case in the explicit example 
in the next section, we can consider examples where A does commute with itself at different times, in which case the 
time-ordered exponential simplifies to an ordinary exponential exactly. Finally, we can compute the Green's function 
numerically and use it in a semi-analytic calculation (remember that the Green's function has to be computed only 
once). That will be worked out in a future publication, though we give some results in section IY Dl 



V. EXPLICIT SOLUTION FOR TWO-FIELD SLOW-ROLL CASE 



In this section we provide an analytic solution for the bispectrum in two-field slow-roll inflation. We assume slow roll 
as in section lTvl but in order to obtain explicit solutions in sections IV Al IV Bl and IV 01 we make the further assumption 
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that all slow-roll parameters are constant. The semi-analytic results of section IV Dl are not bound by this further 
assumption and all slow-roll parameters are calculated numerically for a quadratic model. No slow-roll parameter 
takes values greater than unity. However, a semi-analytic approach is feasible for the most general non-slow-roll case 
and will be the study of a future publication [(|. 



A. Power spectrum 



We now restrict ourselves to just two fields. Moreover, we assume the background values of H and all slow-roll 
parameters, including perpendicular ones, to be completely constant in time whenever they are the leading-order 
(in slow roll) non-zero terms in our expressions. Then we can actually solve the system explicitly, i.e. do the time 
integrals. We start from the equation of motion (|17|l for Vi together with the definitions (J2HJ), or rather from (|19J) for 

the mth order v^™ 1 ^ in an expansion in perturbation orders. We assume everywhere that x > 0. At first order in the 
perturbations this reads as 



iP+A^v^=b^, 



lim vP = 0. 



(43) 



The matrix contains just background quantities, which by assumption are constant. Hence we circumvent the 
issues of non-commutativity and time-ordered exponentials, and we can write down the solution immediately as 



At'e Amt 'bf\t\x) 



(44) 



(In the terminology of scction lTTTl the Green's function is G{t, t') = exp[— A^(t — t')].) The exponential can be worked 
out using its eigenvalues and eigenvectors. When multiplied with its inverse (at a different time) and expanded to 
first order in slow roll, we obtain 



(45) 

(i \ (i-{y/ y 'f) ^(i-(i + f)(yAyr + f(y/y') 3 - x ) % (i - {y/y'f (i + f ) ({y/y') x - {y/y'f- x ))\ 
o {y/y'f 2^((y/yr-(y/y'f- x ) ^ ({y/y'f + f {y/y'Y - (l + | ){y/y'f- x ) 

(i + | )W) x_| W )3- x + (y/y'f-x) 



\o o -x {{y/y'f - {y/y'f - x ) -%{y/y') x + (i + f){y/y'f- x 



J 



Obviously, it is the identity matrix if y' = y. For calculational simplicity here we have defined the new time variable 
y, as well as the relative momentum p, 



II ; 



V2 



fc*i?, c 
V2 ~ V2 ' 



-At, 



- A 



kR 



(46) 



with At* = t—t*, the time since horizon crossing of a reference mode k* , and we have used the fact that k* exp(— 1« ) = 1 
by definition (Atk is defined similarly for the mode k). The fixed reference mode fc* is most conveniently chosen to 
be one of the observable modes, say the one that crossed the horizon 50 e-folds before the end of inflation. From 1231) 
and H38fl and the relation a = kcj (pyH\/2) we find that b\ can be written to leading order in slow roll as 



H 



2V2- 



d 3 k 



1 



2p 2 y 2 e~ p ' y ik.e 1 



ikx 



( ai(fc) 
2^ ± a 2 (k) 
a 2 (k) 
\-Xa 2 (fe), 



c.c. 



( 27r )3/2 k 3/2 

Changing to y as integration variable, we can then do the integral in (|44|) explicitly to find the solution 

d 3 k 



(47) 



k H 



1 



2V2ViJ (2^) 3 / 2 fc 3 / 2 



ih e l * 



l ai+2 n-a 2 ~2^a 2 0\ ( 

2^ ± a 2 

a 2 

V - X a 2 0/ V~ x y 3 ~ x ri 



2 2 

e -P y 



pVr(-i,pV) 
p*y*r(i-f,pV> 



|,p 2 y 2 )/ 



H 



2y/2- 



d 3 k 



1 



(2tt)3/2 fc3/l 



iki e ik x B(k)u(py) + c.c. 



c.c. 



(48) 
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where we have omitted the explicit k dependence of the a's and the final expression defines the matrix B(k) and 
vector u(py). 

It is interesting to look at the time behaviour of l|48(l in more detail. Note that in our leading-order slow-roll 
approximation, differences At in the time variable defined in J2J are equal to differences in the number of e-folds. A 
few e-folds 4 after horizon crossing the vector u(py) in l|48|l can be approximated by (1, 0,p x y x T(l — x/2),0) T . The 
third entry can be approximated even further as just 1 — xAifc, where Atk is the number of e-folds after horizon 
crossing of the mode k 7 and the expression is valid for xA£& sufficiently smaller than unity, but Atk ~ 3. 5 With this 
approximation the solution l|48|l can be written as 



vP(t,x) 



K H 

2V2 Vl 



d 3 k 



( 27r )3/2 fc 3/l 



■ ifc,- e 1 ' 



/ a x + 2fj ± At k a 2 \ 
2f, ± {l- X At k )a 2 
(1 - xAtfe)a 2 
\-x(l - X&t k )a 2 J 



+ c.c. 



(49) 



As expected (see e.g. ^3|) we find that the effectively single-field (ax) component of Q reaches its constant final value 
right after horizon crossing, while the influence of the perpendicular field direction (a 2 , 'isocurvature mode') on Q 
continues to grow with time on super-horizon scales. The velocities Oj and Of are both suppressed by an additional 



slow- roll factor compared to the Q's. In the limit of At k — > oo (or py 
longer valid, the exact result (|48() leads to the limit 



0), where the above approximation is no 



vf\x) 



K H 

2V2 VI 



d 3 fe 



ax 



( 27r )3/2 k 3/: 



■ ih e* 



ifc- a: 



V 



2--a 2 






(50) 



/ 



Hence the expression does not diverge as t grows, but reaches a well-defined value, which is independent of the 
smoothing parameter c. 

Concentrating now on the adiabatic (ei) component of £ = d~ 2 dXi we find from l|48|) to leading order in slow roll: 



^\t,x) = 



K H 

2V2 VI 



d 3 fc 



1 



(2tt) 3 / 2 k 3 / 2 



-p 2 y 2 (At,)„ n.\ , o 7 ? 



ax(k) + 2-!-g{p,x,AU)a 2 {k) 
X 



c.c, 



where y as a function of At* is given in (|46(l . p — fc/fc*, and we have defined 



p x y x T (l - |, 



2 2 
P V 



1 — p x c 



Xo -xAt, 



(51) 



(52) 



where the approximation is good from a few e-folds after horizon crossing. Hence the two-point correlator is given by 

d 3 fc 1 



8 e 



;c (i,i (t,«K (i,i (t,«o) 

or, equivalently, for the power spectrum: 

;C (1)1 (M)C (1)1 (M) 



(2tt) 3 fc 3 



4 e fc 3 



9 (p,X,&t*) 



ik- (x—x ) 



+ c.c. 



( ^) 5 2 0' \ A/ : ; 



(53) 



(54) 



Here we used (II II) to take the average. Alternatively, we could have used H26|l directly. From a few e-folds after horizon 
crossing, exp(— 2p 2 y 2 ) sa 1 and g(p, x, At*) is given by the final expression in i|52[l. so that the power spectrum is 
independent of the smoothing parameter c. Finally, we can compute the adiabatic spectral index using the expressions 
in ^3]> where the Up e in that paper can be read off from Ij54(l . once the transient horizon-crossing effects have 
disappeared, to be 2(fj J - /x)g{p, Xi At*)e 2 , 



n ad - 1 = -4e - 2ryll - 8^ ^— g{p, X , At 



l-g(p, X ,At t 



(^y 9 2 (p,X,AQ 



(55) 



4 For example, 3 e-folds is good enough if c = 3 and x = 0.05, and this result depends only weakly on the values of c and \- 

5 See the previous footnote. A logarithmic dependence on c has been ignored here. For c = 3 this term is 4 times smaller than 
when Atk = 3, and becomes even less important as At^ grows. 
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Second-order solution 



At second order in the perturbations we expand all quantities in A and bi as explained in (I20|l . using (|35|l . resulting 
in the expressions in i|37|) and (|38|) . Remember that superscripts within parentheses denote the order in perturbation 
theory, while the superscripts without parentheses indicate the component of the vector within the field basis as 
defined in (@J. The resulting equation for has the same structure as (@3J), but with a different source term: 



( \ 


V 3x (1) / 



.(1)2 



lim v. 



(2) 



0. 



(56) 



where b\ is the vector obtained by perturbing H, e, rj^, x, and e\ and 62 inside oi\ and a 2 in b\ given in 147|l . 



Explicitly, the right-hand side of equation 156(1 is to leading order in slow roll given by 

f(2e + ^)a l (k)-fj ± a 2 (k) «i(fc) -r^a^k) - (e + ^)a 2 (k) -a 2 {k)\ 



8 e 



d 3 fed 3 fe' ik, e ikx 



( 27r ) 3 



B(k)u(py) 



x (B{ti)u(qy)e ik ' x + 
n 2 H 2 ff d 3 kd 3 k' ik ie ik 



V 

/ 

2(e- 27jU)rr L 


V -fl>i 

c.c. ) + c.c. 



fj-t-aiify + (2£ + ^)a 2 {k) a 2 (k) (e + jjH)ai(fc) - r^a^fe) ai(k) 



J 



2^ 





-2fj A 


5e4 





-6x- 2(e + 7711)7711 +2(77^)2 


-ip 2 ~ (6- 10e + 277ll +3x)?7 J 





6e- 


-3?7 J 



277 



nil 



? JJ (2^) 3 



2p 2 y 2 e- p y B(k) - 3(0,0, 1,0) B(k)u{py)F [B(tf)u(qy)d 



ik -x 



C.C. + c.c. 



(57) 



where we have defined q = fc'/fc*, as well as the matrices -B(fe) and -F in the last expression (the matrix B and vector 
it were defined in H48|) ~). The entries indicated by an asterisk in the matrix B are not given explicitly here, but can be 

read off from H38fl : they do not contribute to Q 1 and Q 2 to leading order in slow roll, because they are one order 
higher than the corresponding entries in the first and third row, after cancellations in the final result have been taken 



into account. The solution for v\ (t, x) is now given by the same expression (|44l) as v\ (t, x), if one replaces b\ in 
that expression by l|57(l . though actually calculating the time integral to obtain a completely explicit expression is 
clearly more difficult. 

To get all the time-dependent terms together, it is useful to change from the matrix notation used above to a 
component notation, as defined in (|17fl . We define the indices a,b,c,d,e,f running from 1 to 4 to label the components 
in the 4-dimensional 9j, £|, 6 2 } space. Moreover, we rewrite the matrix in (|45|l as 



-A(°h„A(°h' 



= K abc w c (y,y'), 



(58) 



/(l, 0,0,0) ±(1,-1,0,0) ^(l,0,-l-|,f) ^(1,-1,-1-^,1 + ^ 



K 



\ 









(0,1,0,0) 2^(0,0,1,-1) 22_(0,l,§,-l-f) 








(0,0,1 + f,-f) 
x(o, 0,-1,1) 



±(0,0,1 



2\ 



-1 



2X • 



(0,0,-f,! 



w(y,y') 



f) / 

which defines the rank-3 matrix K and the vector w(y,y'). Then the solution for v\ can be written as 



( 1 

{y/y'f 
(y/y') x 
V(y/y') 3 -*' 



k 2 H 2 ff d 3 fed 3 fe' ihe ik * 



O") 3 fclfc 



—K abc (B de (k')e ik '- 
1 2 V 



c.c. 



(59) 



00 dv' 

B bd (k) I dy' 2p 2 y'e-P'y^w c (y 7 y')u e (qy') - 3B 3f (k)F bd \ JL w c (y,y')u e (qy')u f (py') 



c.c. 



(To be precise, c, e, and / are actually indices in a completely different space than the other indices, but it is also 4- 
dimensional.) It is useful to consider the contributions of the two integral terms within the square brackets separately, 
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since they have a different origin (cf. (129(1 1. The hrst term is the variation of the stochastic source, represented in 
(|29|l by the X term, and because of the window function the integral only picks up a contribution around horizon 
crossing (although this contribution is time-dependent even later on, because of the dependence on y, not just y', of 
the Green's function). The second term is the variation of the coefficients in the equation of motion, represented in 
(|29|l by the A term, which is an integrated effect up to the end of inflation, and is not present in single-field inflation. 
The leading-order coefficients in front of the first term are first order in slow roll, while the ones in front of the second 
term are second order 6 , however, this can be more than compensated by the larger integration interval. 

In principle there are 80 different integrals here: 16 from the first term and 64 from the second one. Some of the 
integrals can be done analytically, but most have to be studied numerically. However, of those 64 from the second 
term the only integrals that matter are those that are secular, i.e. continue to grow with time (up to a time of order 
At ~ x -1 ), since these will be, roughly speaking, a slow-roll order of magnitude larger at the end of inflation than the 
other integrals. Although we studied all integrals more carefully, one can easily get an idea of which integrals in the 
second term will be secular by looking at the behaviour of the integrand for y' — > (i.e. t — ► oo): only the components 
with e and / either 1 or 3 are secular, since these are close to y'~ in that limit. (Actually the components with e and 
/ either 2 or 4 are zero in this slow-roll approximation, as can be seen from l|48fl). A slightly more careful analysis 
shows that, roughly speaking, the c = 2 and c = 4 components of those terms will be a factor x smaller (one gets a 
3 _1 instead of a x _1 when integrating). Given that B31 = 0, / cannot be equal to 1, so in the end one expects the 
4 integrals in the second term with c and e both either 1 or 3 and / — 3 to be dominant, and that is confirmed by a 
careful numerical study. We denote the (c, e, /) = (1, 1, 3), (1, 3, 3), (3, 3, 3), and (3, 1, 3) integrals in the second term 
within the square brackets of Q59JI by h(p, q, x, At*), hip, q, X, At*), hip, q, X, At*), and Ii(jp, q, x, At*), respectively. 

Regarding the 16 integrals of the first term the following can be said. Because of the y 3 factor in front of the 
c = 2,4 integrals, which cannot be completely canceled by factors coming from the integral, these terms will become 
negligible after just a few e-folds after horizon crossing. Of the remaining integrals those with e = 2 and e = 4 are 
zero. Hence there are also only 4 distinct integrals here that have to be considered: those with c = 1, 3 and e = 1, 3 
Again, these simple estimates are confirmed by careful numerical study of the integrals. We denote the (c, e) = (1, 1), 
(1,3), (3,3), and (3,1) integrals in the first term within the square brackets of by Ji(p, q, At*), q,Xi At*), 
h(p,q,X, At*), and J 4 (p, q, x, At*), respectively. 

Let us now investigate these 8 integrals. Half of them, viz. ^3, I4, J3, and J4, are zero in the limit of t — > 00, but 
decrease slowly enough with time that they should not be neglected at the end of inflation. Three of the integrals can 
be done analytically: 



Jib, g, At*) = / dy'2p 2 y'e 

Jy(AU) 

J 3 (p,g,X,At*) = qXy*(AU) [ dy' <lp>y' e^^Y (l - §, q 2 y' 2 ) (60) 

= -piTftp 2 + ? 2 ) x/2 y x r (1 - f , b 2 + ?V) + 9 Ve- pV r (i - f , 9 y ) , 



ly(At,) P 2 + Q 2 



•2 



J.ip^xAU) = y x (At*) / dy>2pY~ X e- (p2+ ^' 2 = ^— 2 {p 2 + q 2 ) x/2 V x T (l - f , (p 2 + q 2 )y 2 

Jy(At,) P +q V 2 

where y(Ai*) is given in (|46|l . It is also interesting to look at the behaviour of the integrals in the limits of p — > 
(k — * 0) and q — * (k' — ► 0). For p — > all integrals are zero. For q — > only I\, I4, J\, and J4 are non-zero. The 
integrals J\ and J4 are given above, but in this limit also I\ and I4 can be computed analytically (the expression for 
J4 is only valid from a few e-folds after horizon crossing, i.e. for py <C 1): 



I 1 (p,0,x,At.) = r dyW~ 1+X r(l-|,py 2 ) = -(e-*V-pW(l-|,jV)) = - <?(p,X, At*), 
J y (At,) V2 / x v v 2//x 

7 4 (p,0,x,At*) = pV(At*) r ^r(l-^, P V 2 ) - -p x y x \n(py)T(l-f) « (At* - hxp)p x e^ At ' (61) 

i 9 (Ai.) V V 2 ' V 2/ 

where g(p, x, At*) is defined in l|52|) . Using the results discussed in the text above equation (|49|l . one sees that from a 
few e-folds after horizon crossing both start growing linearly with Atk (= At* — \np), although finally the limit 1/x 
is reached for Ji, while I4 goes to zero. 



There are some entries in the product FtdB^e that are first order in slow roll, but these exactly cancel when 1591 is worked out explicitly, 
so that the non-vanishing leading-order coefficients are second order in slow roll. 
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(c) x = 0.05 (d) x = 0.1 



FIG. 1: The integrals h(p, q, x, At,) (small dashes), hip, q, X, At*) (solid), hip, q, x, At*) (large dashes), and /4(p, q, X, At*) 
(dots) plotted as a function of q/p for At* = 50, both for p = 1 (blue, darker) and for p = exp(10) (green, lighter) (i.e. 50 and 
40 e- folds after horizon crossing of the mode k). The smoothing parameter c = 3, although the dependence on c is negligible. 
The different figures correspond with different values of \, as indicated. Note that h and I4 almost coincide in the first plot. 



For q > the four /-integrals have to be evaluated numerically; the resuls are plotted in figure Q]as a function of 
q/p for At* = 50, for various values of the parameters x an d P- Although we will be using the exact numerical results 
for all integrals when plotting the three-point correlator, one can get an approximation by neglecting the x/2 inside 
the gamma function in uzipy') (see (15911 and <|48l) ). Then the /-integrals can be done analytically, with the results 

hip, q, X , A**) « \v X ip 2 + q 2 r x/2 T (|, ip 2 + q 2 )y 2 ) , hip, q, X , At*) « ~ p*q\p 2 + q 2 )^r ( X , (p 2 + q 2 )y 2 ) , 

hip, q, X , At*) » \p*qHp? + g 2 )-*/ Vr (|, + <zV) , hiP, q, X, At*) « \p x y*T (0, (p 2 + g 2 )y 2 ) , (62) 

so that we can make the following estimates: 

h « At fc (1 - i x At fc + i(xAt fe ) 2 ) , J 2 w At fe (l - x At fe + f (xAt fc ) 2 ) , 

7 3 ~At fc (l-|xAt fc + |(xAt fc ) 2 ), J 4 w At fc (l-xAtfe + i(xAt fc ) 2 ), (63) 

for x 1 ^ Atfc, and 

h~X~\ h~\x-\ h~0, h*0, (64) 

for x _1 ^ Atfc. As a rough approximation, they can be taken independent of q for reasonable ranges, say up to 
q/p ~ 100. For I 2 and ^3 this range has a lower limit as well: 100 -1 ^ g/p ;$ 100; they are zero for q = 0. Note that 
these secular /-integrals typically give a result which is of the order of an inverse slow-roll parameter. 
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FIG. 2: (a) The integrals Ji(p, q, At*) (small dashes) and q, X, At,) for x — 0.05 (solid) and x = 0.2 (large dashes) plotted 
as a function of q/p, a few e-folds after horizon crossing of the mode k when the dependence on Atk = At, — lnp has become 
negligible, (b) The integrals ,h(p, q, Xi At,) (solid) and J<i(p, q, x> At,) (small dashes) plotted as a function of q/p for p = 1 
and At, = 50, for both x — 0.01 (blue, darker) and x = 0.05 (green, lighter). 



For the J-integrals the results are much smaller, since the integration interval is restricted because of the window 
function. The integrals J\ and J 2 become completely independent of Atk from a few e-folds after horizon crossing of 
the mode k. Moreover, J\ is independent of x, while J 2 has a relatively weak dependence on x- On the other hand, 
both depend strongly on q/p. They are plotted in figure |2~(a)| The integrals J3 and J4 depend strongly on both q/p 
and x^ifc- In the limit xAifc "C 1 they beco me equal to J 2 and J%, respectively, while in the opposite limit they both 
go to zero. They are plotted in figure 2(b) For q = we have an exact analytic result; for q = p (i.e. k = k') it is 
sometimes useful to have an analytic approximation: 



Ji 



J 3 m i(l - xAt fc ), 



J 4 « i(l-xAt fe ), 



(65) 



for x 1 > Atk, and 
Ji ■■ 



J-2 



J 3 ~o, 



J 4 ~0, 



(66) 



for < At k . 

Having studied all the integrals, we can now work out (|59|l explicitly. We focus on the a = 1 component of vt a , that 

is Q (the adiabatic component of since that is the quantity we want to compute the three-point correlator of in 

(21 1 

the end. The final result for Q ' in the two-field case, in a leading-order slow-roll approximation (constant slow-roll 
parameters) and valid well after horizon crossing, is 



Cf }1 (^) 



8 e 



d 3 fcd 3 fe' 1 
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with oj\ = x(— c + 2?7" — £ J ~/fj ± ) an d ^2 = (e + 77" + (3e — ?j")x + (?7^) 2 defined as short-hand notation. The 
arguments (p,q,x,At*) of the integrals have been suppressed, but it should of course be kept in mind that that is 
where the time dependence resides in this expression. Note that in the single- field limit, where all terms with a 2 
disappear and t) 1 - = 0, we recover the result of p|. For the three-point correlator we need to know £^ 1 = d~ 2 d l Q 2 , 
which is given by the same expression (|67(l . but with e lk x (ifej e lfc x ai(fc) + c.c.) replaced by 

and the same for ct2(k). 



C. Bispectrum 

As in the single-field case, (C*- 2 -* 1 ) is indeterminate. To remove this ambiguity and also require that perturbations 
have a zero average, we define C m = C™ ~~ (C™)- Expanding ( m — Q 1 ^ m + Q 2 " 1 m and switching over to Fourier space, 
we finally arrive at our end result for the three-point correlator (or rather the bispectrum) of the adiabatic component: 
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Again, this result is valid in the two-field case, in a leading-order slow-roll approximation (constant slow-roll parame- 
ters) and valid from a sufficient number of e-folds after horizon crossing that transient effects have disappeared. The 
function g{p, x, At*) is given in (|52|) . %, ipi, and ip 2 are defined in section llVl and u>i and c^2 arc defined below l|67|l . 
Remember that all the integrals and the function g(p,x> At*) depend on the momenta via p — k/k* and q = k'/k* 
and hence are affected by the interchange of k and fc'. In the single-field limit only the first term on the first line of 
(|70|l remains, which agrees exactly with [5j. 

In the limit A: 3 <C fci, k 2 (and hence k\ = — = fc, while we also fix k* — k so that we do not need to write a 
subscript on At), all the integrals can be performed analytically and the result is (leaving aside the overall factor of 
(27r) 3 <5 3 (£ 5 fc s )): 



( 2 ) = «_j_j^ 

8 k 3 kl e 2 I V X 



(2?+r/H) 1 + 4 



(1 



-*At\ 



(71) 



•4( 2_ 

X 



^ (1 - (1 + xAt)e^ At ) + ^ (1 - e~* At ) 

X X 



where the term on the first line within the curly brackets comes from the J-integrals, and the term on the second 
line from the /-integrals. Again, this agrees with the single-field result in the limit fj — > 0. Unlike the single-field 
case, the multiple-field result cannot be expressed in terms of the scalar spectral index and the power spectrum only 
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FIG. 3: The bispectrum 17111 divided by the square of the power spectrum I54H (with two different momenta) in the limit where 
one of the momenta is much smaller than the other two, plotted as a function of fj ± and % f° r At = 50, e = —77" = 0.05, and 
(V2lAc)Viii/(3# 2 ) = -£ x = -(\/2I/ac)VW(3# 2 ) = 0.003. 

(see (|54[1 and (|55|) . and 17] for expressions for the isocurvature and mixing components). Instead of the three-point 
correlator itself, it is actually more useful to look at the ratio of the bispectrum to the square of the power spectrum, 
since that ratio is related to observables like the /nl parameter (more about that later). Dividing H71|) by the square 
of 154fl (one with momentum k and the other with £3) and taking the limit of fj- 1 /x 3> 1, we get for the two opposite 
limits of %At that 



Now if we assume that fj is larger than the other slow-roll parameters, the dominating term in both expressions 
will be the 4(r)- L ) 2 in ipi l|36|) . so that the two expressions in i|72|) will go to 4(77 _L ) 2 At and 8(77 ) 2 /x, respectively. 
Hence, while this ratio of the bispectrum to the square of the power spectrum is first order in slow roll by naive power 
counting (counting 1/At as a slow- roll parameter), as in the single-field case, it can be much larger for models with 
a relatively small x an( i relatively large fj- 1 -. For example, fj = 0.07 and At = 1/x — 50 would already give a ratio 
of more than unity, so that a value about 100 times larger than in the single-field case seems well within range for 
multiple-field models. This is confirmed by the full plot of (|71|l divided by the square of l|54|l as a function of fj 1 - and 
X given in figure [3] It is also interesting to see that in the cases where non-Gaussianity is large, this is caused by 
the /-integrals (i.e. the super- horizon integrated background effects that are absent in single-field inflation): roughly 
speaking it boils down to eJi versus (fj ± ) 2 Ii, which gives e versus the smaller of (?)- L ) 2 At and (fj ± ) 2 /x, either of which 
can easily be two orders of magnitude larger. 

In the opposite limit of k\ = k2 = &3 = k (where we again set fc* = k so that At is unambiguous) we do not have 
an exact analytic result for all integrals, but we can use the approximations i|63l) (|66[1 . We find the following results 
for the bispectrum divided by the square of the power spectrum, in the limit of ry^/x 3> 1: 



If we assume once again that fj is larger than the other slow-roll parameters, the dominating term in both expressions 
will again be the 4(?7 J -) 2 in ipi, so that the two expressions will go to 2(fj ± ) 2 At and 6(ry- L ) 2 /x, respectively. Finally 
we should check that all the limits that produce large non-Gaussianity do not produce an unacceptably large spectral 
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FIG. 4: (a) The bispectrum I69L 170H without the overall (27r) 3 <5 3 (^ s fc s ) factor and divided by the square of the power 
spectrum 154H as defined in 17511 . plotted as a function of the relative size of the three momenta. The sum of the momenta is 
chosen as (fei + ^2 + ^3) = 3fe» with fe» fixed by choosing Ai» = 50. The values for the parameters are e = — fj" = 0.05, r) ± = 0.2, 
X = 0.01, (V2i/K,)Vm/(3H 2 ) = -g x = ~(V2i/ k)V 2 2i/ (3H 2 ) = (V2l/«)VW(3ii" 2 ) = 0.003 (as well as c = 3, although the 
dependence on c is negligible), (b) An explanation of the triangular domain used, defined in 176L with k = ki + ki + kg. 



index at the same time. Fortunately that is not the case: from l|55(l we derive, under the same limits as in (|72J) and 
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After having discussed the various momentum limits, we finally show the full dependence on the relative magnitude 
of the momenta of the bispectrum divided by the square of the power spectrum in figure 4(a) where we did not 
use any analytic approximations for the integrals. To be precise, it is actually the bispectrum given in (|69fl and 
(|70|l . without the overall (2tt) 3 8 3 (J2 s k s ) factor (but taking into account the relation between the momenta that the 
S- function implies), divided by the sum of products of power spectra (|54|l with different momenta, as follows: 



/nl 



(75) 



This quantity can be seen as a momentum-dependent version of the /nl parameter often used in the literature (see 
e -g- Q )- 7 We now choose to be the mode that crossed the horizon 50 e-folds before the end of inflation (i.e. we 
set At* = 50). The function /nl depends on the three scalars fci,A;2, A3, but we can plot it in a two-dimensional 
triangular domain if we fix their sum, which we do by setting {k\ + + kz)/k„ = 3. This convenient way of plotting 
the three-point correlator in a triangle, clearly demonstrating its symmetries, was introduced in 5], and is illustrated 
in figure 4(b) 8 The quantities on the axes are 



k-2 - &3 
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0= -V3 



(76) 



At the vertices of the triangle one of the three momenta is equal to zero. Lines of constant k s are parallel to the 
sides of the triangle (a different side for each s = 1, 2, 3) and k s increases linearly perpendicular to these. At the side 



There is a difference of a factor of order unity between /nl and /nl even in the equal momentum limit, caused partly by the difference 
between f and the gravitational potential $ which was used in the original definition. 
Note, however, that in |^ a different normalisation factor was used. 
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itself the corresponding momentum is equal to half the total sum, (k% + k-i + k%)/2. In the centre of the triangle all 
momenta have equal length. 

The plot in figure 4(a) has been made for all first-order slow-roll parameters equal to 0.05, except rj 1 - = 0.2 and 
X = 0.01, and all second-order slow-roll parameters equal to 0.003. We see that there is a dependence on the relative 
magnitude of the momenta. Though not visible in the figure, this dependence is strongest very near the vertices of the 
triangle, which is the limit of (|71|l . where for this specific example the value 9.4 is reached. Of course logarithmically 
the region near the vertices covers an infinite range of magnitudes in momentum ratios. (The fact that the result is 
largest in the squeezed momentum limit agrees with the findings of [18(-) The value at the centre is 3.7. Assuming that 
a naive extrapolation of this result at the end of inflation to the time of recombination is allowed, so that the quantity 
plotted is indeed comparable to the observable /nl> we see that this model does produce sufficient non-Gaussianity 
to be detectable with the Planck satellite. To compare this plot with the one for the single-field case in [j| one should 
keep in mind that there an additional factor of (2e + ?y") was left out (and there are some differences in the momentum 
normalisation factor, but that does not change the magnitude much), so that the multiple-field result is indeed about 
two orders of magnitude larger. 



D. Comparison with semi-analytic calculation using quadratic potential 

Of course it may be argued that the approximate model considered here is not very realistic, with all slow-roll 
parameters constant with time (in particular rj- 1 and x). We should also stress that the calculation here was made 
under the assumption that \ > 0, which is not true for all models. While we will study more realistic models in great 
detail in a future publication both semi-analytically and purely numerically without any approximations, for direct 
comparison here we present a bispectrum calculation using the Green's function formalism outlined in section ITTT1 We 
have investigated a simple two-field model with a quadratic potential V = ^m\(j)\ + \m%4^ with m\ = 1 • 10 _5 k _1 
(the overall mass magnitude can be freely adjusted to fix the amplitude of the power spectrum). The analytic solution 
(|30H is used as the linear source term in (|23[) . the super-horizon Green's function is then calculated from (1221) and 
(|33|l . and the bispectrum computed from these using (|29|l and 1)3 7|) . We find that for a mass ratio rri2 / ttv\ — 9 and 
initial conditions (f>i = <t>2 = 13k -1 we get relatively large non-Gaussianity: with all momenta equal, that is, at the 
centre of the triangle, the ratio of the bispectrum to the square of the power spectrum is, in the slow-roll limit, 

/nl = 1.8, (77) 

where we have taken horizon crossing to be 58 e-folds before the end of inflation. The ratio of the contribution from 
the /-integrals to that of the J- integrals is 74. This confirms our assertion that the integrated secular terms (the A 
term in $29\ ) subsequent to horizon crossing dominate the contributions to the bispectrum. We note also that the 
spectral index in this model is 0.93, which is observationally acceptable. While the investigation of the quadratic 
model is preliminary at this stage, it is clear that large non-Gaussianity (/nl greater than unity) can be obtained in 
a real multiple-field inflation model. 

Even though the slow-roll parameters are definitely not constant in the quadratic model, we find that the expressions 
in the previous subsection can be used to make a useful approximation of the numerical result. In this case, for 
example, to estimate fj' L we use a representative value and adjust At to reflect the region of support where its value 
is significant. To compare to the quadratic potential result, here we have taken (r)' L ) 2 = 0.73 (its maximum value) 
and At — 1.0 (full width at half maximum), and use the limit 2(fj ± ) 2 At given below J73J). From this we estimate the 
value 1.5, which is relatively close to the numerical value. This seems to indicate that the constant slow-roll, analytic 
results of the previous subsection 9 , while in principle unrealistic, can be used to get a first estimate of the amount of 
non-Gaussianity even in real models with varying slow-roll parameters. 10 



Except expression 1741 for the spectral index, which is a poor approximation in this case. 

Note added: After our investigation of this example in an earlier version of this paper, the non-Gaussianity produced by a quadratic 
potential was also considered by the authors of |19| . They conclude that Jml <S 1 using the '5N formalism' and extrapolating results 
from a potential with (almost) equal masses. They have not computed the general case of unequal masses and assume deviations from the 
(almost) equal mass case to be small. In the case of almost equal masses, which is effectively single-field, we agree that non-Gaussianity 
is small. However, already in ll.'l it was quantitatively shown that even at linear order additional leading-order effects arise in the case 
of unequal masses from the effective coupling between the fields caused by the bending of the trajectory in field space. Moreover, in the 
model we consider here, the dominant non -Ga ussianity is caused by a relatively large fj 1 - , so that a naive slow-roll order counting is not 
valid, a situation where the calculation of llll as we understand it is not applicable. 
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E. Discussion 



Before we conclude, a couple of points regarding the consistency of our approach need to be discussed. The first is 
an inherent limitation of our method in capturing all possible sources of non-Gaussianity, since, by using the linear 
perturbation solutions to source our non-linear equations, we are implicitly neglecting all non-linear interactions up to 
horizon crossing. We believe that this accounts for the small discrepancy between the momentum dependence of our 
single-field three-point correlator p| and that obtained from the tree- level action calculations of when k\ ~ k 2 ~ k 3 , 
whereas in the limit k\,k 2 3> k% the two correlators agree exactly. We surmise that the super-horizon non-linear effects 
described by our method and the horizon-crossing effects we are missing are of comparable magnitude for single-field 
inflation in the equal momentum limit. 

For multiple-field inflation, however, the situation is very different. We can see this by using the quantitative 
results above to interpret our key integral expression for the three-point correlator Ij29(l . In the case where multiple- 
field effects are large (as indicated by the behaviour of fj' L ) it is the perturbation of the long- wavelength evolution term 

in the integrand of l|29[l (represented by Aw c and absent in the single-field case) that dominates over perturbations 

of the stochastic source term which contains the linear perturbations (represented by xim C ; the term that would, 
in principle, be influenced by these horizon-crossing effects). For example, in the case considered in figure 4(a) the 
contribution of the xim C term, which is given by the J-integrals in (|7U|I . is 62 times smaller than the contribution 

from the A^ term at the vertices of the triangle, and 177 times smaller at the centre. But would the xim C term 
be similarly enhanced when taking into account non-linear effects at horizon crossing in the multiple-field case? This 
seems unlikely, given that horizon crossing is only a short transition, while the large effects of the other term are caused 
by a buildup over a significant time interval. Moreover, this question appears to have been answered definitively in 
the negative by recent work j2jj. Generalising for multiple-field inflation, though only up to horizon crossing, it 
shows that these extra contributions remain of the order of small slow-roll parameters, just as in the single-field case. 
In that sense, the papers 0, l20| are important null results which clarify that our approach focusing on non-linear 
super-horizon effects will indeed capture the main non-Gaussian contributions from multiple-field inflation models. 

The second point regards the possible influence of loop corrections to the stochastic picture for generating and evolv- 
ing inflationary perturbations. It is generally accepted within the cosmological community that quantum fluctuations 
can be considered classical for modes which have crossed the horizon, and we explicitly make such an assumption here 
by using classical random fields to set up initial conditions for long wavelengths via the source terms in (|l(Jfl . The 
long-wavelength evolution is then followed by using the classical equations of motion. A natural question to be asked 
is whether loop corrections might play a role in the super-horizon evolution. Recently, the question was addressed in 
21] for a single inflaton field plus a number of non-interacting massless scalar fields. A theorem was proved about 
the growth of loop effects and it was shown that for the theories mentioned loop effects were determined at horizon 
crossing and were subdominant. Since the conditions of the theorem imply fj = 0, these results are not directly 
applicable to the kind of models considered in this paper. However, even if loop effects were to grow with time in such 
models, they would still need to dominate over the classical growth that these models can exhibit in order to interfere 
with the classical picture for the evolution of the perturbations we have developed here. Nevertheless, a definitive 
answer to such matters requires further investigation. 



VI. CONCLUSIONS 

In this paper we have investigated non-Gaussianity in multiple-field inflation using the formalism of 0, Q , em- 
phasizing analytic calculations. That formalism is based on fully non-linear equations for long wavelengths, with 
stochastic source terms taking into account the short-wavelength quantum fluctuations. For analytic calculations 
an expansion of the relevant equations in perturbation orders is necessary. However, it is much easier to derive the 
perturbed equations at second order directly from the non-linear equation of motion for Q than from perturbing the 
original Einstein equations. Of course, in a fully numerical investigation no expansion in perturbation orders has to 
be made; this will be explored in future work. 

We derived two main results in this paper. The first is the general solution for the bispectrum, (|27|l with l|28l) or 129fl . 
Even though this is an integral expression, it will be relatively simple to evaluate in a semi-analytic calculation and 
it yields the full momentum dependence. To achieve this one only needs solutions for the homogeneous background 
quantities in the inflation model, for the linear perturbation variable Q around horizon crossing, and for the 
homogeneous Green's function, as well as expressions for the spatial derivatives of the various coefficient functions. 
The latter can all be computed analytically from the constraint equations l|14|) (|16fl ; for the general two-field case all 
relevant expressions were given explicitly in |3( and this paper. Computing the bispectrum is then just a question 
of performing a few time integrals. An accurate semi-analytic treatment will be the subject of a forthcoming paper 



20 



0, though we do provide the results of a slow-roll calculation for a quadratic potential here. In the present paper, 
however, we have emphasized an example where we could proceed purely analytically. 

In the second part of the paper we studied two-field slow-roll inflation, with the strong leading-order approximation 
that all slow-roll parameters are constant. In this case we could work out the bispectrum explicitly analytically (apart 
from a few integrals that had to be done numerically, although we found analytic approximations in certain limits), 
which is the other main result of this paper, equation H70(l . We found that in this two-field case the bispectrum can 
easily be two orders of magnitude larger than in the single-field case, due to the continued buildup of non-Gaussianity 
on super-horizon scales caused by the influence of the isocurvature mode on the adiabatic perturbation. We note 
that even though the presence of isocurvature perturbations during inflation is crucial, it is not mandatory that they 
survive afterwards. In fact they feed into the adiabatic perturbation and can disappear by the end of inflation (as in 
the cases studied here). On the other hand, if any isocurvature modes do persist at the end of inflation, their fate 
will depend on the details of reheating and further evolution. 

The bispectrum divided by the square of the power spectrum, which can be seen as a momentum-dependent 
generalisation of the /nl observable, can be /nl ~ 0(1)-0(1O), or even larger in extreme cases. If a straightforward 
extrapolation of this result at the end of inflation to the time of recombination is justified, a subject which still needs 
to be studied in more detail, this means that the Planck satellite, and to a lesser extent even the WMAP satellite, 
should be able to confirm or rule out certain classes of multiple-field inflation models. Finally we want to stress 
that beyond estimating the amplitude of the bispectrum, we also give its explicit momentum dependence. While 
this dependence is rather flat for momenta of comparable size, there is a significant difference between more extreme 
momentum limits. 

We believe this paper is a significant step towards providing quantitative and testable predictions of non-Gaussianity 
from multiple-field inflation. Nevertheless there are still a number of issues that remain to be investigated in more 
detail, and on which we are working for future publications. In the first place, we will apply our general solution 
for the bispectrum to more realistic inflation models, particularly those strongly motivated by fundamental theory. 
This will require a semi-analytic treatment, but because we are dealing with integral equations, we believe that the 
strong slow-roll approximations presented here will actually provide reasonable analytic estimates of the exact results. 
As a first step we presented here the results of the semi-analytic slow-roll treatment of an explicit two-field model 
with a quadratic potential. This will be investigated in more detail in 0, lll| , but the results confirm the fact that 
non-Gaussianity can be large in multiple-field inflation models, and that our analytic approximations provide a good 
estimate. Next, it is of course important to study the further evolution of non-Gaussianities after inflation through 
recombination to the present day (see 0,H2 for some work in this direction). In this paper we restricted ourselves to 
computing only the bispectrum of the adiabatic component of even though we have the solution for all components. 
In future work we will investigate isocurvature and mixed bispectra as well. Finally, we will test our results with 
a purely numerical implementation of our formalism, which can also be applied to non-slow-roll models where our 
analytic approximations fail. In this case, the real-space realisations for £ that result allow for other measures of 
non-Gaussianity to be determined, not just the three-point (or higher) correlator. After the disappointing results 
for single-field inflation, primordial non-Gaussianity is now back as an important quantitative tool for confirming 
or ruling out multiple-field inflation models, offering an exciting new window on the early universe as observations 
continue to improve over the next 5-10 years. 
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